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In the framework of minimal flavor violation (MFV), we discuss the decay properties of a su- 
persymmetric scalar top (stop) in the presence of a light gravitino. Given a small mass difference 
between the lighter stop and lightest neutralino and an otherwise sufficiently decoupled spectrum, 
the stop may be long-lived and thus can provide support to MFV at hadron colliders. For a bino- 
like lightest neutralino, we apply bounds from searches in the 77-^^^ channel (ATLAS with lfb~^ 
and D0 with 6.3 fb"^) and give a 5 fb ^ projection for the ATLAS search. 



I. INTRODUCTION 

Supersymmetry (SUSY) |T| is an attractive exten- 
sion of the standard model (SM) . However the simplest 
version of a supersynimetric SM. the minimal super- 
symmetric standard model (MSSM) [2 , does not pre- 
dict a specific flavor structure; all superrenormalizable 
soft supersymmetry breaking terms allowed by gauge 
and Lorentz symmetry as well as R parity are present 
in its Lagrange density |2] . However, it is clear that a 
supersymmetric extension of the standard model must 
have a non-generic flavor structure to be compatible 
with experimental results [3J H] . 

The way the standard model flavor structure is ex- 
tended to the MSSM is not unique, yet a widely dis- 
cussed flavor scheme is minimal flavor violation [5] 
(MFV). In MFV the standard model Yukawa cou- 
plings are promoted to spurion fields transforming un- 
der the SM fiavor group to restore the SM's flavor 
symmetry. If all additional flavor structure of a new 
physics model can be understood as higher dimen- 
sional flavor invariant operators including these spuri- 
ous and the model's flelds, the model is called MFV. 

As the LHC is running and eventually will flnd su- 
persymmetry, it will be challenging to investigate the 
flavor structure at a hadron collider due to the detec- 
tors' limited flavor identiflcation abilities and the com- 
plexity of the recorded events. In [B] it was pointed out 
that the third generation's squarks decouple from the 
flrst two generations in MFV. As a result, a light stop 
can be long-lived decaying through the flavor changing 
neutral current channel [T, "5^ 

ii ^ cxI, (1) 

if all flavor diagonal channels are kinematically closed. 
{ti denotes the light stop, Xi the lightest neutralino 
and c a charm quark.) An observation of long-living 
light stops thus would hint in the direction of MFV. 
In MFV, the coupling Y between ti, c and Xi is 

Y cx XlV^bVtl (2) 



where A;, and Vij denotes the bottom Yukawa coupling 
and elements of the Cabibbo Kobayashi Maskawa 
(CKM) matrix respectively. The precise value of Y 
depends on the stop left-right composition, the neu- 
tralino decomposition, and on a numeric factor stem- 
ming from the MFV expansion; see Ref. |2] for details. 

In t9j it is shown that the average transverse im- 
pact parameters for the stop decay products can be 
expected to be 0(1800 /im) for stop lifetimes of the or- 
der of ten ps in the production channel pp — > tttiti [W\ 
. When both top quarks in this channel decay leptoni- 
cally, the pair of same signed leptons in the flnal state 
allows to separate the signal process from its SM back- 
ground; however, the small leptonic branching ratio of 
top quarks suppresses this process so that, after ap- 
plying all kinematic cuts, only few events are left in 
this channel. Ref. jTT] proposes an alternative collider 
signature assuming stop pair production in association 
with one hard jet. Demanding a minimum transverse 
momentum of 1 TeV for the additional jet, the whole 
parameter region consistent with electroweak baryoge- 
nesis can be probed. Ref. |T2] considers an analogous 
process, stop pair production in association with two 
b-jets. However, these two studies do not consider 
stops in the MFV framework. 

If we consider local SUSY instead of a global im- 
plementation of SUSY, we can have distinct collider 
signal signatures with little SM background: In lo- 
cal SUSY, a massive gravitino emerges in the super- 
symmetric mass spectrum [13]. Its interactions with 
other particles are severely suppressed by the reduced 
Planck mass 

mpi = (87rG'Ar)-5 2.4 x lO^^ GeV, (3) 

where Gn is Newton's constant. Depending on the 
exact breaking mechanism in the hidden sector, the 
gravitino can be very light. A light gravitino inter- 
acts through its goldstino components with couplings 
proportional to 

(rn3/2mpi)~\ (4) 



2 



where TO3/2 is the gravitino mass. 

In models of gauge mediation [14t418|. TO3/2 is gener- 
ally much smaller than the sparticle mass scale; thus, 
the gravitino is the lightest supersymmetric particle 
(LSP). Its goldstino interactions are enhanced and can 
be of the electroweak order. 

Consequently, the lightest neutralino decays via 

X? ^ XG, (5) 

where X denotes a photon, Z, or a Higgs boson; G 
denotes a gravitino [TO] . If X is a photon (7) , this de- 
cay leads to very clear collider signatures with high pT 
isolated photons plus missing transverse energy (^t) 
stemming from gravitinos leaving the detector unseen. 

Several studies with light gravitinos at hadron col- 
liders were performed in the past. Ref. [19 21j con- 
sider the diphoton plus channel at hadron col- 
liders. In Ref. |22| [23] the authors examine a stau 
NLSP and a gravitino LSP. A sneutrino NLSP and a 
gravitino LSP scenario is investigated in Ref. |24l E5] . 
Ref. IH] considers the discovery potential of a neu- 
tralino NLSP and a gravitino LSP at the Tevatron, 
where they consider a general decomposition of the 
lightest neutralino. In Ref. [19^ the authors consider a 
light stop NNLSP and a light neutralino NLSP and a 
gravitino LSP. Ref. [2T investigate a stop NLSP and a 
gravitino LSP scenario for the Tevatron as well as the 
LHC. A chargino NLSP and a gravitino LSP is consid- 
ered in Ref. [28] . Depending on the size of the grav- 
itino mass, non-pointing photons can be measured. 
The discovery potential of sparticle decays with a fi- 
nite decay length are investigated in [251 - 1^ . A re- 
cent experimental search for sparticles with finite de- 
cay lengths is published in [32'. 

In this paper we investigate the parameter region 
where the decay in Eq. (jlj is dominant in the context 
of a light gravitino, and how the stop masses are con- 
strained in this framework by recent collider searches, 
assuming that ti, Xi, and G are the only light su- 
persymmetric particles. We discuss the stop and Xi 
decay patterns in section |ll] In section |III| we apply 
collider bounds in the 77^7- channel [551 The 
cuts adapted from the experimental studies are dis- 
cussed in the two appendices. 



II. DECAY PATTERNS 
A. Light stop decays 

Here we discuss possible decay patterns of light 
stops and their implications on the parameter space. 
We start with decays of the light stop via Yukawa and 
gauge couplings and then discuss direct stop decays 
into a gravitino. 

The lighter stop mass eigenstate ii is the lightest 
squark state in many supersymmetry breaking scenar- 
ios. On one hand, the large top Yukawa coupling can 



induce a sizable left-right mixing in the stop sector 
leading to light ti masses. On the other hand, if the 
soft squark mass terms are unified at a high scale, 
the stop mass terms are prominently reduced by the 
top Yukawa coupling in the running of the renormal- 
ization group equations to the electroweak scale [35] . 
In this case, ii is mostly a SU{2) singlet state since 
its mass term does not receive any contributions from 
SU{2) gaugino loops. In addition, right-handed stop 
loop contributions to the rho parameter are sufficiently 
suppressed [39 . 

Surprisingly light masses of 0(100 GeV) are still 
consistent with experimental searches for stops decay- 
ing to cXi at the Tevatron if the mass splitting 

Am = ruf^ - m^o (6) 

is smaller than « 30 GeV [3D;. 

Since we want the fiavor changing decay in Eq. Q to 
be the dominant decay and the stop to be long-lived, 
we must ensure that the potentially dominant decays 
ii — >• txi, ii — >• bxi, and ii — > bx^W are kinematically 
closed or sufficiently suppressed, 

nif^ < TOfc + m^o + mw, < mt + m^± . (7) 

Depending on the chargino and slepton 
masses (> m^^), the four-body decay ti — > bivx\ 
may be dominant if Am exceeds a few 10 GeV. For 
small stop neutralino mass splittings, the tree-level 
four-body decay is strongly phase space suppressed 
and Eq. remains the dominant decay mode [3]. 

The MFV decay width can be written as 

in the limit of Am ^ m^^ [6J. Thus, if this decay is 
dominant, the stop lifetime is governed by AmK. The 
kinetic energy of the hadronic stop decay remnants 
depends on the size of the mass splitting Am. In a 
study with same-sign leptons [Sj , a major reduction of 
event numbers due to a minimal pt cut on these low 
energy decay products has been found. 

In the presence of a light gravitino, restrictions on 
Am and Y are not sufficient to guarantee the domi- 
nance of the FCNC (fiavor changing neutral current) 
decay mode Eq. ([ij. If AmY is too small, the two- 
body decays of the stop into the gravitino, 

ti ^ Gc, h -> Gt, (9) 

and — if the resonant decay to tops is kinematically 
closed — the three-body decay mode of the stop, 

fi ^ GbW, (10) 

can have a significant contribution to the full stop de- 
cay width. These decays can invalidate our assump- 
tion, that the stop dominantly decays via a FCNC 
process with a finite impact parameter. 
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For the flavor diagonal two and three-body decay 
channels, the decay rates are |19| |4T| 



ti lifetime [ps] 



r{ii tG) = 

T{ii W+bG) = 



1 



(1 



2\4 



(Ua) 



3847r2 sin^ ( 



2 2 
?T1 x3 1 ^2 



W '"-PI 



• [ICLP + |C^|V(X2^,X2)] , (lib) 

where the gravitino mass is neglected in the phase 
space integrals. (Xerm and dy/ denote the fine-structure 
constant and the Weinberg angle, respectively. Fur- 
ther x-w = "tn-w l"fntx = TTit/Tn^^ with the top 
mass mt- cl and cr parametrize the and tu con- 
tribution to ti] due to the 3rd generation's decou- 
pling in MFV the other squarks' admixture is small, 
i.e. jcip -t- |c_Rp « 1. The functions I{xy^,x1) and 



are phase space integrals and can be found 



in |3T]. Note that Eq. (lib) does not comprise the 



finite top width and diverges at the top mass thresh- 
old. We use this formula for m^^ < mt only. As the 
three-body decay proceeds trough a virtual top quark, 
its rate is largest for a right-handed stop, because the 
chirality fiipping top mass dominates the propagator. 

For a bino-like Xi, the flavor structure of the ti — 
G — c coupling stemming from the MFV expansion is 
the same as in the ti — — c coupling, thus the decay 
rate for ti — >■ Gc can be written as [9^ 



T{ii Gc) 



Y' 



ti 



487r 'Tipjm3^2 ' 



(12) 



where Y' and Y are related by a factor dependent on 
the stop composition as Y comprises the hypercharges 
of the left- and right-handed stop flelds. The factor is 



Y 




(right-handed ii) 
(left-handed ti) 



(13) 



with the SM U{1) coupling g' and the left-handed 
(right-handed) stop hypercharge Yq — g (|). 

We show the branching ratio B{ti — >■ ex?) and the 
stop lifetime in the m^-Y plane for three different 
masses of a right-handed stop in Fig.[l]using the decay 
rates ([8|, (11), and (12). To generate the plots, we 
keep Am fixed at 10 GeV and use sin^ 9w = 0.23, 
aem = 1/128, and mt = 173 GeV. The plot for left- 
handed stops does not differ significantly from the one 
shown. 

Due to the m^^ dependence of the decay widths 
in Eqs ([9]), ( |To| ) and the weaker dependence of 

T{ti — > ex?), the minimal gravitino mass necessary to 
account for a sizable B{ti ex?) increases with larger 
stop masses. The m^^ dependence of r(ii ex?) also 
causes the smallness of the shifts of the lifetime regions 
in Fig. [T] to larger Y values when the stop mass is in- 
creased. 
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FIG. 1: MFV light stop branching ratios (left, blue) 
and stop lifetimes (right, pink), respectively for rai_^ — 
150 GeV, 300 GeV and 500 GeV (from top to bottom), in 
dependence on the gravitino mass rriQ and the MFV cou- 
pling Y. In the branching ratio plots, the shaded areas 
correspond to B{ii cx?) > 0.01, 0.1, 0.9 (from light to 
dark hue), whereas, in the lifetime plots, the shaded areas 
correspond to rj^ > 0.01, 0.1, 1, 10, 100 ps (from light to 
dark hue). We chose Am = 10 GeV. 



As it is clearly visible from Fig. [T] very small val- 
ues of Y < 0(10"^) and at least gravitino masses of 
0(0.1 — IkeV) are required in addition to the mass 
hierarchy 



(14) 



for the stop to be long-lived and to decay dominantly 



to ex?. 



As the charmed hadron produced in the decay has 
a macroscopic lifetime of O(lps) itself, however, a 
macroscopic stop lifetime might turn out to be acces- 
sible experimentally only if it exceeds this timescale. 



B. NLSP neutralino composition and decays 

Neutralinos are mass eigenstates of the C/(l) gauge 
fermion (bino), the neutral SU{2) gauge fermion 
(wino), and the neutral up- and down-type Higgs 
fermion (higgsino). The fields' individual contribu- 
tions to Xi as well as the neutralino mass spectrum 
depends on the bino mass Mi, the wino mass M2, 
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the Higgs mixing parameter /i and the ratio between 
the up-type and down-type Higgs vacuum expectation 
value (VEV) tan /3. If x" is the NLSP and G the LSP, 
decays via x\ — > XG, where X is either the pho- 
ton, the Z boson, or a neutral Higgs boson. Branch- 
ing ratios into the various decay channels are fixed by 
phase space suppression factors and the decomposi- 
tion of the lightest neutralino. General formulae for 
the decay widths are given in Refs [191 ES] . 

In the previous subsection, we argued that the mass 
of the light chargino x J , a mass eigenstate of charged 
winos and higgsinos, has to be larger than the light 
stop mass in order to suppress the flavor diagonal stop 
decay to Xib. This requirement cannot be satisfied if 
the Xi is wino-like, i.e. if M2 ^ Mi, as in this case 
both the Xi ^^^d the x^ mass are « M2. The mass 
splitting between wino-like xf and Xi is of the order 

of [42], given Mi ^ and thus is extremely 
suppressed for > few 100 GeV. 

Similarly, if Xi is higgsino-like, Xi ™d Xi have 
masses of the same order of magnitude given by /i. 

The mass splitting Am^±^(, is of the order of ^ 

for 1^1 < Mi,M2. 

If Xi is bino-like, its mass is Mi approximately, 
while the xf mass is given by \fi\ or Af2- As the mass 
gap depends on two different supersymmetric mass pa- 
rameters, it can be sizable depending on the details of 
high scale physics. 

While in a mass region close to the Z mass, also 
a higgsino-like Xi ma-Y respect the anticipated mass 
hierarchy in Eq. ( 14 1 , we focus on a bino-like Xi in 
discussing experimental bounds as 1) in the bino case 
the mass hierarchy can exist over a large stop mass 
scale and 2) binos have a large branching fraction to 
photons. For m-^^o > mz and negligible phase space 

suppression, the branching ratio is B{xi — ^ 7G) ~ 
cos^ 0w ■ This value is obvious as G is a gauge singlet 
and 7 a mixed state of the hypercharge gauge boson 
and the neutral SU (2) gauge boson where the mixture 
is parametrized by the Weinberg angle. Including the 
phase space suppression from mz and assuming that 
the higgsino sector is decoupled, the bino decay rates 
are [Ulig 



nxi -> 7G) - cos^ 

r(x; ^ ZG) = 



'W 



m-0 

Xi 

'iSnmy^mpY 



sm Vy/- 



'3/2"'Pl 



Xi , 



(15a) 



(15b) 



For reference we show the bino-neutralino lifetime 
as a function of the lightest neutralino mass for grav- 
itino masses 1, 10, 100, 1000 eV in Fig. |2] 




100 150 200 250 300 350 400 450 500 550 
m^o [GeV] 

FIG. 2: Lightest neutralino lifetime as a function of the 
neutrahno and gravitino mass. 

III. COLLIDER BOUNDS 

When ti, Xi, and G are the only light supersymmet- 
ric particles, stops are dominantly produced in pairs, 
both at pp and pp colliders. 



pp titl, 



PP 



titl 



(16) 



The production cross sections are given in Fig. [3] as a 
function of the stop mass for the LHC at 7 TeV as well 
as for Tevatron and are calculated with Prospino | 44| 
using the built-in CTEQ6.6M parton distribution 
functions (PDFs). We also show the next -to-leading 
order uncertainty by varying the factorization scale 
(fip) and the renormalization scale {fiR) between ■^iti^_^ 
and 2m while keeping 11 ^ equal to /jp. 

Given a bino-like x?, the final state signatures of a 
decay chain via Eq. ([ij) and Eq. (|5| are 



77ccGG, jZccGG, ZZccGG. 



(17) 



In general, with a small mass gap Am, the charm jets 
are too soft to be useful for event selection. Thus con- 
strains on our parameter space can stem from searches 
for an excess in the 77^^, ^Z^rp and ZZ^j^ chan- 
nels. As binos dominantly decay to 7G, the SM back- 
ground is negligible for energetic photons, and large 
and high px photons axe efficiently identified in 
multipurpose detectors, we concentrate on the jj^x 
channel in this work. Several experimental searches for 
the diphoton and channel have been published |33f - 
[57] . So far, no excess above the SM expectation has 
been found. 

In the following, we present exclusion limits in 
the stop-gravitino mass plane derived from the latest 
ATLAS search in the j^^x channel for a luminosity 
(£) of 1.07 fb~^ [31]. We derive also bounds from the 
D0 search with £ = 6.3 fb~^ [33], and give a £ = 5 fb~^ 
projection for the ATLAS bound. 

The dominant SM background with originating 
from the hard process stems from W + -f, W + jets. 
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in bin ri from 
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400 500 
mr [GeV] 



600 700 



FIG. 3: NLO total cross section for titl production at the 
Tevatron and at the LHC with i/s = 7TeV. The colored 
bands correspond to scaling the unified scale /ip = I-Ir = 
by ^ and 2. 



and W/Zj{e)^{e) production. Here, electrons/jets are 
misidentified as photons. SM background with 
from mismeasurements emerges from multijet produc- 
tion and direct photon production. 

A supersymmetric background can only arise from 
XiXi P^'ir production as all other colored sparticles, 
sleptons, heavier neutralinos and charginos are as- 
sumed to be heavy and thus will have a negligible 
contribution. However, also the XiXi production cross 
section is severely suppressed. In the limiting case of 
vanishing higgsino admixture to Xi, the cross section 
vanishes even at ©(al^a^). Consequently we do not 
take XiXi P^ir production into account in the follow- 
ing. 



A. Calculation of exclusion limits 

To constrain the stop mass, the gravitino mass, and 
the MFV coupling Y, we calculate cr^^j. , the total cross 

section for tit\ production, in a grid of the light stop 
mass rrij^ using Prospino. Note here that the light 
stop mass is the dominant SUSY parameter in the 
cross section [JJ, both for pp and pp initial states. 

For each stop mass in the grid, we generate 
100 000 iiil pair events with pythia 6.4.25 [46] us- 
ing the CTEQ6.6M [15] parton distribution functions. 
With the hadron level events we simulate the effi- 
ciency/acceptance for the 77ccGG' final state in the 
ATLAS and D0 analyses employing a slightly modified 
version of Delphes 1.9 (l^ISS]- The photon energy, 
and therefore our signal's detection efficiency, depends 
on the mass splitting Am. As in the previous sections, 
we fix Am = 10 GeV. In appendices |X] and |B] we 
give details on the cuts adopted from the experimen- 
tal studies and on further simulation parameters. As a 
result of this simulation step, we obtain efficiencies e„ 
in bins of and can calculate a signal cross section 



e„6(ti ^ X?cr6(x? ^ G7) aij,. (18) 



Using Eqs (15 1 to calculate S(xi — >■ G7), we fi- 
nally employ the CL^ method |48l 115] |57j to calculate 
the 95% exclusion limits for B{ti — >■ Xic)- Those are 
depicted in Fig. [4| where we use the measurements 
plus background predictions of the experimental stud- 
ies as enlisted in Tab. [l] When calculating the exclu- 
sion limit, we treat the errors on the luminosity and 
the background as Gaussian nuisance parameters, see 
Tab. |T| but do not take into account theory uncertain- 
ties stemming from scale variations and the choice of 
PDF sets. 

The projection for the ATLAS study with a lumi- 
nosity of 5fb^^ is calculated using the prescription 
of I50i. 





Bin [GeV) 


Observed 
events 


SM bgd 
events 


D0 [33 


35-50 


18 


11.9 ±2.0 


(6.3±0.4)fb-^ 


50-75 


3 


5.0 ±0.9 




> 75 


1 


1.9 ±0.4 


ATLAS [34] 


> 125 


5 


4.1 ±0.6 


(1.07 ± 0.04) fb~^ 









TABLE I: ATLAS and D0 measurements and background 
(bgd) predictions. 



B. Numerical analysis and discussion 

As can be seen in Fig. |4] the ATLAS search (solid 
pink curve) gives a bound on B(ti — >■ Xi^) for stop 
masses up to 560 GeV. For a luminosity of 5fb^^, 
this mass is projected to 660 GeV. For larger masses, 
a dominant tree level FGNG decay ti — >■ XiC is in 
agreement with the measurement. As can be seen from 
Fig. m the 1.07 fb-i ATLAS data already excludes a 
larger parameter region than the 6.3 fb~^ D0 data |58] . 

For each stop mass, the bound in Fig. [1] can be 
mapped to a bound on the maximal gravitino mass 
for given values of Y using Eqs ^, and p). 

We plot these bounds for Y = 10"'^, 10-^ lO'^n 
Fig. [sj As Eqs (111 do not provide the correct stop 



width for masses in the threshold region close to the 
top mass, we exclude this region from the mapping 
and interpolate our result in the region mt ± 30 GeV 
(mt = 173 GeV). 

The mass difference Am enters the bounds in Fig. |5] 
through Eq. ^ and through the hardness of the pho- 
ton pt spectrum. For small changes in Am, the latter 
dependency can be neglected, and the bounds depend 
on the product (FAm) only. Therefore the bounds for 
other viable values of Am can be derived from those 
shown for Am = 10 GeV by rescaling Y. 



6 




\oO 200 300 400 500 600 700 



mi, [GeV] 

FIG. 4: Maximal branching ratio for i\ — > XiC in depen- 
dence of the stop mass given a bino-like xi and a fixed stop 
neutrahno mass difference of Am = 10 GeV. Blue curve 
and blue shaded area: D0 study P^. Pink curves: ATLAS 
search [34], where the solid curve corresponds to the mea- 
surement, and the dashed line corresponds to a £ = 5fb~^ 
projection. The areas above the curves are excluded at 
95% CL within our simplified analysis. 

Obviously, the smaller we choose Y the larger 
can be to generate a branching ratio below a certain 
value. The bound on 7713/2 resulting from the mapping 
in Fig. [5] varies over a wide mass range and thus poten- 
tially implies different regimes of <i and x? lifetimes; 
therefore, we show slopes of fixed values: For the Xi 
lifetime, we show slopes for Ifs, Ips, and Ins (black 

5/2 

dotted) following approximately 7713/2 oc m-^ . For ti, 
we show slopes for 1/5 ps, Ips, and 5ps (black solid). 
At Y = 10^^, the stop lifetime is below Ips already 
induced by ti — ^ j^c irrespective of the gravitational 
decay channels; therefore, only the 1/5 ps slope can 
be drawn here. Similarly, the stop lifetime drops be- 
low 1/5 ps for stop masses smaller than w 230 GeV 
for Y = 10^^, as T{ti — ex?) oc ?7i^^ for fixed A?7i. 
For both smaller values of F, the stop lifetime slopes 
shown only depend weakly on Y because the corre- 
sponding total stop widths are dominated by the grav- 
itational decay ti — ?> tG resp. ti — ^ W~^bG. 

Fig. [5] shows that the gravitino mass region pro- 
moted in [3] as the region where the decay ti — > XiC 
dominates for stop masses between 100 and 150 GeV 
is now disfavored. More generally. Fig. [5] allows to 
discuss two regimes of different orders of ii and x? 
lifetimes: 

• In the stop mass region where B{ii — Xic) is 
bounded, i.e. for m^^ < 500 GeV, the grav- 




[GeV] 

FIG. 5: Pink curves: The maximal gravitino mass respect- 
ing the ATLAS measurements |34| in dependence of the 
light stop's mass. The solid curve corresponds to the mea- 
surement for C = lfb~^, and the dashed line corresponds 
to a C — 5fb~^ projection. Blue curve: C = 6.3 fb~^ 
D0 search |33| . Gravitino masses above the curves are ex- 
cluded at 95% CL. The black dotted curves are slopes of 
fixed Xi lifetime (Ifs, Ips, and Ins from bottom to top) 
for 771^0 = mj^ — 10 GeV. Along the solid black lines, the 
stop lifetime is 0.2 ps, Ips, and 5ps. 
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100 200 300 400 500 600 700 
[GeV] 

FIG. 6: White area: the available ratio of xi lifetime over 
ii lifetime for Y = 10"®, Am = 10 GeV. Blue area: ex- 
cluded by the D0 data [33]. Pink area: excluded by the 
1.07 fb"^ ATLAS data [HI]. Pink dashed line: 5fb"^ pro- 
jection for ATLAS. Grey area: invalid. 



itino channels have a significant contribution to 
the stop decay. If this contribution is dominant, 
both — the Xi ^^nd the ti decay — are governed by 
the same coupling ~ l/rriy^- As the stop decay 
width is suppressed by phase space {tG channel) 
or top propagator {GW~^b channel), the stop life- 
time is expected to be larger than, or at the same 
order of magnitude as, the neutralino lifetime in 
this region. 



where B{ii — > Xic) is al- 
for mi^ > 500 GeV, the 



In the mass region 
lowed to dominate, i.e. 
phase space suppression of the stop decay width 

Eq, 



(11a 



is less pronounced; thus, the stop and 
neutralino gravitational partial decay widths are 
of the same order of magnitude. Consequently, if 
— > cxi is the dominant stop decay in this mass 
region, the stop lifetime is significantly smaller 
than the neutralino lifetime. 

The relation between the stop and neutralino life- 
times (tj^ and r^o) described above is summarized in 
Fig.|6]where we plot the allowed ratio of both lifetimes 
within the bounds. Along the solid black line bound- 
ing the grey area, the contribution of Eq. ([8| to the 
stop decay width vanishes; thus, this line represents 
the smallest value the ratio can have. The plot is gen- 
erated for Y — 10^^; however, only the leftmost parts 
of the exclusion areas depend weakly on the specific 
value of Y . For smaller masses, the bounds are driven 



by the phase space dependence of Eqs and ([15| . 

Note that for large ii or x? lifetimes, care must be 
taken in the interpretation of the bounds in Figs |5] 
and [6j The selection criteria for photon candidates in 
the ATLAS publication are chosen for prompt pho- 
tons |5T1 15^ . Also, more explicitly, the D0 study 
requires that photon candidates point to a recon- 
structed primary vertex. We assumed in our calcu- 



lation that all signal photons fulfill these criteria as 
if they were prompt. In a more realistic simulation, 
for increasing neutralino lifetimes, the signal photons' 
selection efficiency should decrease. For longitudinal 
neutralino decay lengths of 0(1000 mm) ATLAS sim- 
ulations show that for several photon selection crite- 
ria jSn] used in [M] the efficiency drops from 0(85%) 
to 0(55%) |53]. Consequently, for lifetimes > ©(ns), 
we overestimate the number of photons accepted, and 
the bounds on presented should be regarded with 
care in this lifetime region. 

While for a bino-like neutralino, the 77^7- channel 
offers the highest sensitivity for setting mass limits, it 
may be difficult to measure the neutralino lifetime in 
this channel due to lack of photon tracks. To construct 
the photons' impact parameters, CMS focuses on con- 
verted photons in a 2.1 fb~^ search in the 77^7- +jets 
channel in [54j. As pointed out in [31 , the x° ~^ ZG 
channel may be used to investigate the neutralinos' 
lifetimes, as the Z's decay products allow to recon- 
struct the Xi's trajectory. 



IV. SUMMARY 

In SUSY models with MFV, the third generation 
of squarks decouples from the other two generations. 
This decoupling opens the opportunity to support the 
MFV hypotheses with the measurement of a macro- 
scopic stop decay length if the stop decay can only 
proceed through a generation-changing channel due 
to kinematic constraints |B]. In this work, we investi- 
gate the decay of a stop into a charm and a bino-like 
lightest neutralino with a subsequent bino decay to a 
photon and a light gravitino. 



X?H 7G)c, 



(19) 



given a sufficiently small mass difference between ti 
and 

While a macroscopic stop decay length serves as 
a hint for MFV, the neutralino's decay to a photon 
leaves a distinct signature in LHC and Tevatron de- 
tectors offering a good signal isolation. 

Assuming that the remainder of the SUSY spectrum 
is decoupled, we find that the ATLAS search in the 
"i^$T channel based on a luminosity of 1 fb~^ jM] im- 
plies a bound on B{ti — ^ ex?) for stop masses up to 
560 GeV, see Fig.ji] In a 5fb^^ projection, the bound 
is raised to 660 GeV. 

We find that stops lighter than ^ 400 GeV are still 
compatible with the 77-^T searches; however, in this 
region, a significant fraction of the stops decays into 
gravitinos and quarks of the third generation. Here 
the stops are expected to have larger lifetimes than 
the lightest neutralinos; though, a macroscopic stop 
decay length is governed by the gravitino mass scale 
TO3/2 and is no hint for a decoupled stop fiavor mixing 
structure. 
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Stops heavier than ^ 500 GeV can dominantly de- 
cay through the decay chain in Eq. (19 1 and eventu- 
ally support MFV if they are a long-lived. As the 
lifetime depends on both the stops's flavor mixing and 
the gravitino mass, it can serve as a hint to MFV only 
if the contribution of the gravitino modes to the total 
decay width is negligible. This is assured if the neu- 
tralino lifetime is much larger than the stop lifetime. 

We have discussed a split spectrum where the only 
light particles are a light gravitino, a light bino, and a 
light right-handed stop. It may be possible, however, 
to separate supersymmetric background processes in a 
scenario with additional light sparticles by vetoing on 
additional jets and/or isolated leptons. 
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Appendix A: D0 cuts 

Irrespective of the very clear signal event structure, 
we employ the fast detector simulation Delphes as a 
framework to calculate the signal efficiency. For D0 we 
use a simplified calorimeter layout composed from cells 
of dimension O.lx |^ in rjx (p space covering \rj\ < 4.2 
and (j) G [0, 2tt[ where rj denotes pseudorapidity and (jj 
the azimuthal angle. 

Jets are constructed employing the iterative mid- 
point algorithm with R = 0.5. 

We adopt the cuts of [33J by requiring 

• at least two isolated photons with > 25 GeV 
and |?7| < 1.1, 

• the azimuthal angle between $x and the hardest 
jet, if existent, is < 2.5, 

• the azimuthal angles between $!rp and both pho- 
tons are > 0.2, 



• > 35 GeV. 

Photons must have 95% of their energy deposited in 
the electromagnetic calorimeter. For a photon to be 
isolated, the calorimetric isolation variable / defined 
in [33] must fulfill / < 0.1 and the scalar sum of 
transverse momenta of the tracks in a distance of 
0.05 < R < 0.4 from th e photon m ust be smaller 
than 2 GeV. Here is R = ^ Ac/)^ -f- Ary^^ where A(/) is a 
track's azimuthal distance from the photon and Ary is 
the corresponding distance in pseudorapidity. 

Note that nearly all signal events fulfill the isolation 
criteria hinting that the hadronic stop decay products 
are well separated from the photon stemming from the 
subsequent neutralino decay. This is expected as the 
neutralino decay products 7 and G are massless and 
thus the photon can have a large relative to the 
stop fiight direction. 

Also note that we assume that a primary vertex 
can be reconstructed, and that the photon trajecto- 
ries point to this vertex. The latter assumption should 
be regarded with care when the neutralino lifetime is 
large. 



Appendix B: ATLAS cuts 

We use the default ATLAS detector layout imple- 
mented in Delphes and apply the following simplified 
cuts: 

• At least two isolated photons exist with pT > 
25 GeV and \t]\ < 1.81, but outside the transition 
region 1.37 < I77I < 1.52, 

• $T> 125 GeV 

The ATLAS study employs a tight photon selection 
criterion on photon candidates, where the efficiency to 
identify a true (prompt) photon is approximately 85% 
in the kinematic region considered |51l 152] . We mimic 
this selection criterion by removing photons from our 
Monte Carlo sample with a probability of 15%. 

We consider a photon to be isolated if in a cone of 
width R = 0.2, the scalar Et sum is less than 5 GeV. 
Here we sum the Et of the Delphes calorimeter cells 
and exclude the cell the photon is mapped to. 

As for the D0 case, nearly all signal photons fulfill 
the isolation requirement. For larger Xi masses, the 
main reduction of signal event numbers stems from the 
85% photon selection efficiency. 
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